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Abstract. We would like to investigate the information contained in our 
observations and to what extent each of them contributes individually to con- 
straining the physical parameters of the system we are investigating. To do this, 
we present a study involving the technique of Singular Value Decomposition us- 
ing as a simple example a detached eclipsing binary system. We intend to apply 
an extension of this technique to asteroseismic measurements of Delta Scuti stars 
that are members of eclipsing binary systems. 



1. Introduction 

Asteroseismology is the study of the oscillations of the stars. It is the only 
method which allows us to see beyond the photosphere of a star, and learn about 
their interior conditions. We are interested in 5 Scuti stars, whose pulsation 
amplitudes are large enough to be detected from the Earth. However, they 
present a physical problem: not all of the oscillation modes that theory predicts 
are excited nor observed, which leads to various interpretations of the frequency 
spectra and thus a poorly-constrained solution (mass, age etc.). If we study 5 
Scuti stars in eclipsing binary systems we could accurately derive their stellar 
masses and radii. Then we could constrain some of the parameters of the system. 
This in turn would lead to an easier identification of the modes of the stars, 
allowing us to better determine the real parameters of the system. Combining 
information from these different sources has been suggeste d by many authors, 
and s ome studies have been conducted with this in mind ( Aerts fc Harmaned 



I2nn4l : lEodnguez et al.ll2nn3 'l. In fact, they agree that this is the only way to 
advance if we would like to test stellar structure and e volution theor i es. 

Following a technique similar to that conducted bv lBrown et al. I (119941 ) ■ we 



would like to investigate how the combination of information from the study of 
binary systems with oscillation frequencies constrains the solution parameters. 
In particular we would like to understand by how much each observable (radial 
velocities, oscillation modes etc.) influences each of the parameters of the system 
(mass, age etc). To approac h this problem w e use the technique of Singular 
Value Decomposition (SVD) ()Press et al.lll986h and apply it to a simple case of 
a detached eclipsing binary system. 
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Table 1. Initial parameters and observational errors 



Parameters 




Observables 


0"i 


RAiRcs) 


0.850 


KA{kms-^) 


0.90 


i?s(R©) 


0.095 


A'B(kms-i) 


0.90 


in 


89.300 


Dp (mag) 


0.02 


T2/T1 


0.300 


LF(hrs) 


0.20 


a(AU) 


0.060 


LT(hrs) 


0.20 




0.970 


P(yrs) 


5x10^ 


Mb{Mq) 


0.100 


Ds/Dp 


0.012 



2. Mathematical Background 



For a more detailed discussion we refer the reader to iBrown et all ^99^). Let 
Bi,i = 1,M be the observables of some system, and Pj,j = 1,N with M > N 
be the variables (parameters) that define this system. There exists some model 
or list of equations that when given Pj, will result in Bi. In real situations we 
know Bi and we would like to extract Pj. If the system or model were linear, 
this would be achieved relatively easy by solving a system of M linear equations 
in N variables. However, normally this is not the case, and so by linearizing the 
system around a reference set of parameters, Pjo, we can write 

^ SB 
j=i J 

where BiQ are the observables resulting from the initial parameter estimates Pjo, 
6Pj = Pj — Pjo and each of the derivatives, is evaluated at Pj = PjQ. The 

parameter fitting problem involves finding 6Pj . Given the real observables f3i and 
after some substituting our problem is mathematically equivalent to minimizing 

X 
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\D6P - (55||2 where SBi = (A - Bio)/ai, and 



I dBi 

D = —^ (2 

(Ji dPj ^ ' 

The solution is 5P = VW^^U'^dB, where D = UWV^ is the singular value 
decomposition of D. It is precisely this decomposition that provides us with 
information regarding the origin of the errors in the parameter solutions. In the 
absence of the real observables Pi we can still construct the derivative matrix D 
to study the transformation between the observable and the parameter space. 



3. Application to a Detached Eclipsing Binary System 

Given photometric and radial velocity measurements of a detached eclipsing 
spectroscopic binary system, we would like to obtain the following fundamental 
parameters: masses of both components and Mb, radii of both stars Ra and 
Rb, effective temperature ratio Tb/Ta, inclination of system i and separation of 
both components a. In order to obtain these parameters we need the following 
observables: Ka and Kb are the semi-amplitudes of the radial velocity curve. 



Relationship Between Observations and Stellar Parameters 
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Figure 1. Information from the decomposition matrices of the derivative 
matrix D. Left: Three parameter solutions showing the extent to which the 
parameters are constrained. Center: Three observable vectors quantifying the 
amount of information contributed by each observable. Right: The effect of 
increasing the precision of some observations by a factor of two on each of the 
solution parameter errors. 



From the photometric light curve we obtain the primary eclipse depth Dp and 
with the secondary Ds, a depth ratio D^/Dp. We also observe a total duration 
of eclipse Lt, and a duration of eclipse minimum, Lp (see assumptions below). 
The photometric period P is obtained from long term observations. Let us also 
assume for simplicity that there is no eccentricity e = 0, that neither of the 
components have stellar spots, that the secondary eclipse is total and that both 
stars are well isolated such that there is no mass transfer between them. 

The model that relates these given observables to the parameters is a set of 
equations that can be easily derived. After obtaining observations with errors cjj, 
we make some initial parameter estimates Pjq (Table 1). We can then calculate 
the matrix D (equation l^j) . Using SVD, we decompose D to obtain its singular 
values and its corresponding orthogonal matrices. Figure ^ shows the first three 
solution vectors (left panel) and observable vectors (center panel) of D. The 
vector corresponding to the largest singular value (smallest error) is the top 
solution vector. According to the figure, the best-constrained parameter is the 
separation between both components (left panel, top row vector). There is 
also a small effect from this singular value on both masses. If we inspect the 
observable vector (center panel, top row vector), we can see clearly that P is 
the most reliable observation (corresponding to the largest singular value). We 
know that P is directly related to a and Ma + Mb through Kepler's Law giving 
a well-determined a. In this case Ma and Mb are mainly constrained by the 
radial velocity measurements (left and center panel, bottom row vector). 

The second solution vector shows that both Lt and Lp are also some of 
the more important observables in constraining the solution. Figure^left panel 
shows that both Ra and Rb correspond to this singular value, implying that Lp 
and Lt are mainly responsible for the reduction in the errors of both radii. In 
fact, Lp and Lp contribute also to i but to a lesser extent (sixth singular value). 

Suppose now that the errors in our measurements (cij) were to decrease 
by a factor of 2. By decreasing each of the cij individually, we may investigate 
how each of the errors on the parameters (e^) improves (noting that combining 
all CTj should result in a corresponding reduction of 0.5 in all the e^). Figure 1 
right panel indicates the decrease in each Cj as a fraction of the original error 
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(cjo — where ej comes from the covariance matrix (Press et al. 1986), 

and Cjo represents ej given the original (Xjo, and Cji are the resulting ej given a 
decrease in a particular ai. We show the effects of reducing just a^p (solid) and 
gka (dash). 

The new a^p is solely responsible for the reduction of both e/j^ and eng (a 
total of 0.5). If we were interested in improving the errors on both radii, then 
we know from this that wc need only improve the precision of Dp (note with 
(Tio Lp and Lp were important for Ra and Rb)- Similarly, ctka responsible 
for all of the reduction of cmb- However, (Tk^ is also partially responsible for 
reducing both ^"^^ (by 0.2 and 0.1 respectively). This indicates that if 
we were interested in improving the precision of Ma , we would need not only an 
improvement of (TKa-, but a combination of improved cTj (in this case (JKbi ^Lp 
and ctlj,). We should also take into account the actual values of the errors, for 
example, if e^j^ improves from 0.004 to 0.002 Mq (which is this case), we could 
also conclude that these observations may not be worthwhile. So we must look 
at these results cautiously, and we must also emphasize the sensitivity of these 
results to uiq. 

If we could decrease each Ui by 3, 4, 5... times, we could determine at 
what error level this observable becomes important in constraining a particular 
parameter, given all the other (jj. This could indicate if it were worthwhile 
obtaining better observations. It may also help determine the precision of our 
observations necessary to obtain a solution to the system with errors to within 
a certain %. Studying the transformation between parameter and observable 
space may allow us to reach these conclusions. 

4. Conclusions 

What we have discussed is just a simple example to illustrate what we would like 
to achieve in the case of a more complicated system, such as a detached eclipsing 
binary system where one of the components is a 5 Scuti star. Unfortunately, 
for that case we can not write simple equations to relate the observables to the 
parameters, and for this reason a study such as this one is beneficial. 

We saw briefly the capacity of SVD to try to explain the relationship be- 
tween the errors in each of the observables to the errors in each of the param- 
eters. We saw how the effect of an increase in the precision of the observables 
can impact in a rather complex form on the resulting errors of the parameters. 

Ideally we would like to observe a real eclipsing system with a 5 Scuti 
component. From studying this transformation we hope to be able to optimise 
our observations to obtain the best constrained parameters. Doing this would 
allow us to solve this system to the required precision in order to test theories 
of stellar structure. 
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